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CN ■ Abstract 



We consider two particular ID quantum many-body systems with local interactions related 
to the root system Cat. Both models describe identical particles moving on the half-line with 
non-trivial boundary conditions at the origin, and they are in many ways complementary to 
. each other. We discuss the Bethe Ansatz solution for the first model where the interaction 

potentials are delta-functions, and we find that this provides an exact solution not only in 
the boson case but even for the generalized model where the particles are distinguishable. 
In the second model the particles have particular momentum dependent interactions, and we 
find that it is non-trivial and exactly solvable by Bethe Ansatz only in case the particles 
' are fermions. This latter model has a natural physical interpretation as the non-relativistic 

. limit of the massive Thirring model on the half-line. We establish a duality relation between 

nI ' the bosonic delta-interaction model and the fermionic model with local momentum dependent 

interactions. We also elaborate on the physical interpretation of these models. In our discussion 
the Yang-Baxter relations and the Reflection equation play a central role. 
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1 Introduction 

\ Quantum mechanical models with interactions are, in general, very difficult to solve, but there 

exist a few important cases where exact solutions are available, allowing them to be understood 
completely. A prominent example is the delta interaction in one dimension which, in the simplest 
two-particle case, is defined by the Hamiltonian 



X 
J3 



H = -dl + c5{x) (1) 

where c is a real coupling constant and x G M the relative coordinate of the two particles, x = 
xi — X2. This latter model is popular because it allows for an explicit solution by simple means: 
since the delta interaction is restricted to a; = 0, it only manifests itself in the non-trivial boundary 
conditions for eigenfunctions ij){x) of H, 

^(0+) = ^-(-0+) 
^'(0+)-V;'(-0+) = cV(0+), (2) 

and these can be easily accounted for (we write ■i/'(=tO"*') short for the left-and right limits 
lim^-jo ^/'(is;), and similarly for the derivative V')- The natural generalization of this model to 
an arbitrary number N of identical particles defines a prominent exactly solvable quantum many- 
body system which, in the boson case, was solved by Lieb and Liniger [1] and, for the general case 
of distinguishable particles, by Yang [2] in a seminal paper where the Yang-Baxter relations first 
appeared. 
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Interactions localized at points have been studied extensively using the mathematical theory 
of defect indices; see [3] and references therein. From these studies it is well-known that the delta 
interaction is only one of many possible local interactions, and a general such interaction can be 
characterized by four real coupling parameters. This can be easily understood as follows: for a 
ID Hamiltonian H = —d^ + v with an interaction v localized at x = all eigenfunctions il){x) 
should be smooth everywhere except at a; = 0, and {Hip){x) = —'ip"{x) for non-zero x. Requiring 
H to be self-adjoint leads to the following condition, 

/ dx ('^ff{x)ij{x) -^^ij"{x)) =0 (3) 

J\x\>0 ^ ' 
for arbitrary wave functions and or equivalently 

- ^V''].=o+ = [?V' - ^V'']x=-o+ • (4) 
General boundary conditions are of the form 

t/,(o+) = niiV'(-o+) + ui2V''(-o+) 

V'(0+) = U21^(-0+)+^X22V''(-0+) (5) 

(and similarly for 0, of course) and are thus parameterized by four complex parameters Ujk 
which, when imposing Q), are reduced to two complex, or equivalently, four real parameters. The 
boundary conditions in Eq. ((2)) are obviously contained in this class of boundary conditions, but 
there are others, most prominently 

^'(0+) = v'(-o+) 
^^(0+) - V(-0+) = AV''(0+) (6) 

which often has been referred to as delta-prime interaction; see e.g. Section 1.4 in [3]. Recently it 
was shown that these latter boundary conditions arise naturally from the Hamiltonian 

H = -dl + \dj{x)d^ (7) 

where the second term has a physical interpretation as a local interaction depending also on the 
momentum p = —\dx [4] . The A^-body generalization of this model is exactly solvable by Bethe 
Ansatz in the indistinguishable particle case when the particles are either bosons or fermions [4,5] 
but, different from the delta interaction case, not in the general case of distinguishable particles [4]. 
Still, this model is complementary to the model with the delta interactions for at least three 
different reasons [4]: firstly, for indistinguishable particles, the delta interaction model is known 
to be interesting only for bosons (since the delta interaction is trivial on fermion wave functions), 
whereas the j5(5p-interaction is trivial for bosons and non-trivial for fermions. Secondly, while the 
delta interaction model for bosons can be obtained as the non-relativistic limit of the quantum 
sine Gordon model, the p(5p-interaction model for fermions naturally arises as the non-relativistic 
limit of the massive Thirring model. Thirdly, there exists an interesting weak coupling duality 
between the fermionic p(5p-interaction model and the bosonic delta-interaction model. 

As is well-know, exactly solvable many-body systems of particles moving on the full real line 
are naturally associated with the root system A^v-i, and they often allow for extensions to other 
root systems such that the exact solubility is preserved [6] . An early example was given by Gaudin 
who solved the Cm root system variant of the delta interaction model for bosons [7], while the 
general case of this model for arbitrary root systems and distinguishable particles was treated by 
Sutherland [8]. As pointed out by Cherednik [9], models related to the root system Cn describe 
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interacting particles on the half line, and the exact solubility requires the so-called Reflection 
equation to be added to the Yang-Baxter relations. The Reflection equation has played a central 
role in many exactly solvable systems with a boundary; see e.g. [10] and the review [11]. 

In this paper we consider the Bethe Ansatz solution of the Cm versions of the models discussed 
above. Similarly as for the ^tv-i case we find that the delta interaction model is exactly solvable 
in this way even for distinguishable particles, but for the model with momentum independent 
interactions we obtain its exact solution only for indistinguishable particles. We also elaborate on 
the physical interpretation of these models as describing particles on the half-line with non-trivial 
boundary conditions at the origin. 

To be more specific, the models we discuss in the paper are defined by the following Hamilto- 
nians, 

N N 

H = -^dl^+ 2ci ^[5{xj - Xk) + 5{xj + Xk)] + C2 ^ 5{xj) (8) 

j=l j<k j=l 

(delta interactions) and 

N 

j=l j<k 

N 

+{dx, + 9..J5(x, + xk){d,^ + d,,)] + 4A2 ^ d,^6{x,)d,^ (9) 

i=i 

(local momentum dependent interactions). For simplicity we assume all coupling constants posi- 
tive so that there are no bound states. Mathematically, the model in Eq. Q is the Cat variant of 
the model solved by Yang [2] , and Eq. © defines the C n variant of the model discussed in [4] . 

The plan of the rest of this paper is as follows. In Section 2 we consider the Cat delta-interaction 
model, starting by deriving the boundary conditions and thus turning the Schrodinger equation 
HiIj = Eip into a well-defined mathematical problem. We proceed to the Bethe Ansatz solution 
of this model where the Yang-Baxter relations and the Reflection equation play a central role. 
We conclude the section by elaborating on the physical interpretation of this model. In Section 3 
we discuss the Cm variant of the interaction model, in large parts paralleling our discussion 
for the delta-interaction in Section 2. We also present a duality relation between the fermionic 
p(5j>-interaction model and the bosonic delta-interaction model. Appendix A gives some details on 
the verification of the Yang-Baxter relations and the Reflection equation. Appendix B contains a 
few mathematical facts about the Weyl group of Cat, and Appendix C gives some details on the 
physical interpretation of these models. 

2 Delta-interaction 

In this section we provide the exact solution of the Cm delta-interaction in the case of distinguish- 
able particles and elaborate on its physical interpretation. 

2.1 Boundary conditions 

The Hamiltonian Q defining the Cm delta-interaction model is only formal, and to determine its 
eigenfunctions we must first convert it into a set of boundary conditions. 

For completeness we start by discussing the Hamiltonian H in Eq. (^, which can be regarded 
also as the one-particle case of the Hamiltonian in Eq. (jHl), -/V = 1. The first step to find the 
eigenfunction ip of H \s to note that the equation HTp = Eip for all x is equivalent to —ip" = Eip 
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for X ^ together with the boundary conditions in Eq. These boundary conditions are 

obtained by integrating the equation Hip = Eip twice: first from x = —0+ to x > and then once 
more from x = —0~^ to x = 0"*" yields the first condition in Eq. 0, and integrating from x = — 0"^ 
to x = 0"'' yields the second condition in Eq. Thus in this case there are two regions free of 
interactions, x < and x > 0, linked to each other by the boundary condition at x = 0. 

For general N, the interaction-terms of the Hamiltonian H in Eq. © are restricted to xj = itx^ 
and Xj = for 1 < j < k < N , and the eigenfunctions ijj oi H therefore obey the simple equation 

\ 

d"^. + E I ij{xi, . . . , xat) = for Xj / ±Xk and xj / 0, (10) 

and for each of the boundaries of the interaction free regions one gets a pair of boundary conditions 
similarly to the ones for = 1, 



[dxj - dxM\xj=±Xk+o+ - i^x, - dxk)i^\xj=±xt,-o+ = 2ciV'U,=±Xfc-o+ (11a) 

V'Uj=+o+ = i'\xj=-o+ 

dx,il^\xi=o+ - dxjil^\xj=-o+ = C2V'Uj=o+ (lib) 



(these conditions are obtained by a straightforward generalization of the = 1 argument above, 
using d^. ± d^^ = 2dx^±Xk)- 

Obviously there are now many more regions free of interactions. One such region is < xi < 
X2 < . . . < Xat, and all others are obtained from this by permuting the particle labels, j pj 
with p e Sjy (= permutation group), and/or reflecting some of the coordinates, Xj — > —Xj. Thus 
all regions free of interactions can be characterized as follows, 

< CJiXpi < a2Xp2 < . . . < (TNXpN < oo (12) 

where aj = ±1 and p € S^; we will refer to these regions as wedges. It is important to note that 
they can be labeled by elements Q in the group 

Wn := (Z/2Z)^ X Sn (13) 

where the first factor corresponds to the reflections while the second factor corresponds to the 
permutations of the coordinates, 

XQj = cTjXpj for Q = (cji, . . . , (Tn',p) S Wn with Uj G {±1} and p £ Sn- (14) 

In the sequel we will therefore use the following convenient notation for the wedges, 

Aq : < XQi < XQ2 < ... < XQN (15) 

with Q £ Wn. It is interesting to note that the group Wn is isomorphic to the Weyl group of the 
root system Cat; see e.g. [12]. 

2.2 Bethe Ansatz 

Using the boundary conditions deduced in the previous section we now proceed to determine all 
eigenfunctions of the Cn delta-interaction, starting by recalling the physical motivation of the 
Bethe Ansatz below. For that we first consider the Hamiltonian H in Eq. In this case there 
are eigenfunctions ip{x) = exp(iA;x) for x < which are equal to a particular linear combination of 



4 



exp(ifcx) and exp(— for x > 0. This can be interpreted as scattering by the delta interaction 
oc 6{x) where a plane wave is partly transmitted and partly reflected. Regarding H in Eq. 
as a two particle Hamiltonian with x = xi — X2 the relative coordinate and k = {ki — /c2)/2 
the relative momentum, we can interpret this very fact as scattering of a plane wave solution 
exp(iA:iXi + 1/^22^2) into a linear combination of this wave and another one where the particle 
momenta ki and k2 are exchanged, exp(iA;2a;i + 1/01X2)- This suggests that an eigenfunction tp 
of the A^-particle Hamiltonian in Eq. (jHJ which is equal to a plane wave exp(i^^]^ kjXj) in one 
wedge Aq 1)15(1 will be transformed into a linear combination of plane waves exp(i fejXj) in 

any other wedge where kj = CFjkpj, with aj = ±1 resulting from the interactions oc 5{xj) which 
can invert momenta, kj — > —kj^ and p E Sn resulting from the interactions oc 8{xj — x^) which 
can interchange momenta, kj ^ kg. 

We thus see that the group in Eq. ((T^ naturally appears again, kj = kpj for some P £ Wn, and 
the discussion above suggests the following Bethe Ansatz for the eigenfunctions of the Hamiltonian 
H in Eq. ©, 

V'(x) = Ap{Q) e'^'^-^« for < xqi < xq2 < • ■ ■ < xqn (16) 

PgWn 

with X = (xi, . . . , xat) and kp -xq = X^jLi ^Pj^Qji for Q S Wn- The corresponding eigenvalue 

is obviously E = Ylf=i ^j- 

One now has to take into account the boundary conditions in pial b). For each Q G Wn, the 
wedge Aq ((T3|) participates in boundaries: xqi = XQ(j_|_x) for z = 1, 2, . . . (A^ — 1) and xqi = 0, 
and for each of these boundaries we will get two conditions. More specifically, the boundary 
at XQi = XQ(j+i) is between the wedges Aq and Aqj- where Tj € Wn is the transposition 
interchanging i and (i + 1), and the conditions implied by Eq. ()llaj) for j = Qi and k = Q{i + 1) 
are 

Ap{Q) + AptM = Ap{QTi) + ApxXQTi) 
i{kpi - kpi^,+^^)[ApT^{QTi) - Ap{QTi) + AptM - Ap{Q)] = 2ci[Ap{Q) + AptM]- (17a) 

The boundary at xqi = is between the wedges Aq and Aqp^ with Ri G Wn the reflection of 
the first argument, i.e., XR^j = Xj for j / 1 and —Xj for j = I, and the conditions at xqi = 
implied by Eq. ()llb|) for j = Ql are, 

Ap{Q) + Apr, (Q) = Ap{QRi) + Apr, {QRi) 
ikpi[Ap{Q) - ApR,{Q) + Ap{QRi) - ApR,{QRi)] = C2[Ap{QRi) + Apr,{QRi)]. (17b) 

We thus have 2A(2^A!)^ linear, homogeneous equations for the (2^A!)^ coefficients Ap{Q). 
The following beautiful argument due to Yang [2] shows that this system of equations has enough 
non-trivial solutions and, at the same time, gives a recipe to compute all the Ap{Q). 
For that it is important to note that Wn plays a third role: defining 

{R)q,Q' = ^Q',QR (18) 

one can write 

Ap{QR)= {R)q^q,Ap{Q') = {RAp){Q) (19) 

Q'&Wn 

where the first equality is a trivial consequence of the definition, and in the second we interpret 
{R)q^qi as elements of an n x n matrix R with n = 2^ N\ the rank of Wn- These matrices 
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obviously define a representation R ^ R of Wj\f acting on the coefficients Ap{Q). It is worth 
noting that this is identical with the so called (right) regular representation of Wn- 

We can therefore insert Apnp.{QTi) = {TiApT.){Q) in Eq. (|17aj) . and by a simple computation 
show that these latter equations are equivalent to 

Ap = yi(A;p(i+i) - kpi)ApT, (20) 

where we have introduced the operator 

iufi + cj 

Yi{u) = — (21) 

m — c\ 

and interpret Ap as a vector with 2^ N\ elements Ap{Q). In the same way we can rewrite the 
conditions in Eq. H17b|) using ApJi-^^{QRl) = {RiApji-^){Q), 

Ap = Z{2kpi)ApR, (22) 

with the operator 

Z{u) = i!f^l±f^. (23) 
m - C2 

It is well-known that the group Wn is generated by the reflection Ri and the transpositions Tj 
(see e.g. page 21 in [13]). Thus one can use the identities in Eqs. ((JOJ, (|^^ and (|19() to calculate 
recursively all coefficients Ap{Q) from Ai{I) using the operators Z and Yi above. It is important 
to note that there is a possible inconsistency arising from the fact that the representation of 
an element P in Wn as a product of the Tj's and Ri is not unique. However, any two such 
representations can be converted into each other by using the defining relations of the group Wn^ 

TiTi = 1, TiTj = TjTi, for \i - j\ > 1 

TiTi+iTi = Ti+iTiTi+i (24a) 
RiRi = 1, RiTi = TiRi, for i > 1 

RiTiRiTi = TiRiTiRi. (24b) 

Thus no inconsistency can arise provided that 

ApT^T^Q) = Ap{Q), ApnT.iQ) = ApT^T^iQ), for \i - j| > 1 

^pr,T,+iT, (Q) = ApT^^^T^T,^^ (Q) (25a) 
ApR,R,{Q) = Ap{Q), Apr^tAQ) = ApT^R^iQ), for i > 1 

ApR^TiRiTiiQ) = ApT^R^TiBiiQ) (25b) 

for all P,Q£ Wn. Using the recurrence relations H2U|) and (|22() one finds that these conditions 
hold true if and only if the following operator relations are fulfilled, 

Yi{-u)Yi{u) = I, Y^{u)Yj{v) = Yjiv)Yi{u), for K - j| > 1 

Y^{v)Yi+i{u + v)Y^{u) = Y^+^{u)Yi{u + v)Yi+i{v) (26a) 
Z{-u)Z{u) = I, Z{u)Yi{v) = Yi(v)Z{u), for i > 1 
Z{2v)Yi {u + v)Z{2u)Yi {u - v) = Yi {u - v)Z{2u)Yi {u + v)Z{2v) (26b) 

for all real u and v. The validity of this system of equations is necessary and sufficient in order for 
the Bethe Ansatz above to be consistent and the model at hand to be exactly solvable. The first 
three relations are the so called Yang-Baxter relations, and the last one is the Reflection equation. 
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The validity of these relations for arbitrary Tj and Ri can be checked by straightforward but 
somewhat tedious computations (of course, the validity of the Yang-Baxter relation in this case is 
known since a long time [2], and this seems to be the case also for the Reflection equation [8, 10], 
but for completeness we provide the essential steps in the verification in Appendix A.l). 

Thus the Bethe Ansatz H16() is consistent even in the general case of distinguishable particles, 
and we can calculate all coefficients Ap from Aj using the recurrence relation 

Ap = Wp{k)Ai (27) 

where Wp{k) is a product of the operators li(/cp(j+;^) — kpi) and Z{2kpi) obtained by using 
repeatedly ((201) and ((22|) . 

Interesting special cases of this solution are when the particles are indistinguishable, i.e., when 
the particles are fermions of bosons. In the former case Tj = — /, and Eq. ()21() implies Yi{u) = —I 
independent of the coupling constant ci. This shows that the delta interaction is trivial for 
fermions. In the boson case we have Tj = +/, and Yi{u) is a non-trivial phase. As discussed 
in more detail below, there are two different boson cases with different physical interpretations, 
namely Ri = —I and Ri = +1. 



2.3 Physical interpretation 

As is well-known, the Cat delta- interaction model describe interacting particles on the half- line 
with particular boundary conditions at the origin [9]. However, the general solution of the Cat 
delta-interaction model without any restrictions includes many more eigenfunctions than any 
model on the half line, and the relation between these models is therefore not completely obvious. 
In this section we discuss the relation of these models in more detail. We also give a physical 
interpretation of the boundary conditions which occur as limits of particular external potentials 
restricting the particles to the half line. 

As discussed in Appendix B, in any irrep of the group Wn the reflections Rj of the particle 
coordinate xj are represented either by i?j = or — 1. For simplicity we now discuss in more 
detail the cases where all Rj are the same, either +1 or —1, which from a physical point of view 
are the most interesting cases. As discussed in Appendix B, these irreps of Wn can be rather 
easily understood since they are related in a simple way to irreps of Sat. Thus we can impose the 
following restriction on the eigenfunctions ip of the Hamiltonian in Eq. Q, 

{Rjip){xi,. ..,Xj,.. .,xn) = tp{xi, -Xj, . . .,xn) = ±ipixi, ...,Xj,.. .,xn). (28) 

With that assumption we can restrict ourselves to xj > 0, and the boundary conditions in Eq. 
((TTa|) and Eq. (|llb|) reduce to 

V'Uj=Xfe+0+ = '^\xj=Xk-0+ 
(9x,. - 9^J'0U,=a;fc+O+ - - 9^fe)V'U,=Xfe-0+ = 2ciV'U,=a;fc+0+) (29a) 

and 

25^^VU,=o+ = C2V'U,=o+ for Rj = +1 ^^g^^ 
^^1x^=0+ = for Rj = -1 

respectively. These are exactly the boundary conditions obtained from the Hamiltonian 

N 

Ho = -Y.dl+2ciY,Hxj-Xk) (30) 

j=l j<k 
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describing particles on the half-line, Xj > 0, and the boundary conditions at the origin given in 
Eq. (l29b|) . 

It is also interesting to note that these later boundary conditions are obtained by allowing the 
particles to move on the full line, Xj S M, and adding a particular external potential V{xj) to 
the Hamiltonian in Eq. pOf) which effectively constrains the particles to the half line Xj > 0. To 
be specific, these potentials are given by 

y..^iVoQi-x) + {c2/2-VVo)6ix) if% = +l ^ 
\ Vo&{—x) ii Rj = —1 ' 

where Q{—x) is the Heaviside function (equal to one for x < and zero otherwise), and one 
has to take the strong coupling limit Vq — > oo: as shown in Appendix C, in this latter limit the 
eigenfunctions of the Hamiltonian Hq + yixj) on the full line, Xj G M, coincide with the ones 
of Hq on the half-line, xj > 0, and the boundary conditions in Eq. H29b|l . 

As already mentioned, the most important cases in applications are the ones we have considered 
here, i.e., where all the Rj are the same. Nevertheless it would be of interest to consider the 
implications of allowing the Rj to take on different values, in effect dividing the particles into two 
groups distinguished by their interactions with the boundary. 



3 Local momentum-dependent interaction 

In this section we discuss the model with local momentum dependent interactions defined by the 
Hamiltonian in Eq. ®. While most of our discussion is in parallel with the one for the delta 
interaction model in the previous section, we find that the Bethe Ansatz is consistent only for the 
indistinguishable particle case. We also present a duality relation between fermionic variant of 
the model here and the bosonic model discussed in the previous section. 



3.1 Boundary conditions 

We start by considering the = 1 Hamiltonian H in Eq. ((Jj). To obtain the corresponding 
boundary conditions we first integrate from x = —0~^ to x = 0^ which yields the first condition 
in Eq. and integrating from x = —0^ to x > and then once more from x = — 0"*" to x = 0"^ 
yields the second condition. As in the delta interaction case, the eigenfunctions ^ of if are then 
determined by these conditions together with the equation —ip" = Eip for x ^ 0. We note that 
the wave functions ^(x) on which H in Eq. ^ is defined can be discontinuous at x = 0, and to 
make sense of the interactions we have implicitly used a regularization which amounts to replacing 
ip'{0) by [^'(0+) -|- ip' {—0'^)]/2 (this is discussed in more detail in [4]) 

It is straightforward to generalize this argument to the A^-particle case. Similarly as in the delta 
interaction case one finds that the eigenfunctions ip of the Hamiltonian in Eq. ® are determined 
by Eq. (|10)) together with the boundary conditions 

'^\xj=±Xf,+o+ - i^\xj=±Xk-o+ = 2Ai((9^^^ - 5^jVU,=±Xfc-o+ (32a) 
V'U,=o+ - V'U,=-o+ = 4A25^,V'U,=o+- (32b) 
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3.2 Bethe Ansatz 



We now discuss the Bethe Ansatz for the eigenfunctions of the Hamiltonian H defined in Eq. Q. 
Obviously much of what we said for the delta interaction case carries over straightforwardly to the 
present case. Due to the different boundary conditions in Eqs. (|32al b) Eqs. (|17al b) are changed 
to 

■i{kpi - kp(^i+i))[ApT^{QTi) - Ap{QTi)] = \{kpi - A;p(i+i))[Ap(g) - ^pt.(Q)] 
Ap{QTi) + ApT^QT,) - Ap{Q) - AptM = 2\i\{kpi - A;p(i+i))[Ap(Q) - AptM] (33a) 

iA:pi[^p(Q) - ApR,{Q)] = ikpi[ApR,{QRi) - Ap{QRi)] 
Ap{Q) + ApR,{Q) - Ap{QRi) - Apr,{QRi) = A\2ikpi[Ap{QR{) - ApR.iQR^)]. (33b) 

We now also use Eq. (|19() to convert these in to the recurrence relations 

Ap = Y,{kp^^,-kp^)ApT,, (34) 

and similarly 

Ap = Z{2kp,)ApR, (35) 

where now 

iui - l/Xifj 

= —■ TTT — 36 

m — 1/Ai 

and 

mi - 1/A2^i . . 

m - I/A2 

As in the delta interaction case these relations allow to recursively compute all coefficients j4p in 
terms of Aj, and the conditions for the absence of inconsistencies are identical to ()25al b) of the 
delta-interaction case, leading to the Yang-Baxter relations (|26ajl and Reflection equation ()26b|) 
but now with the operators (|36l) and (|37jl . In contrast to the delta-interaction case, we find that 
these consistency relations are valid only if Tj = ±7 for all i (see Appendix A. 2 for details). We 
thus conclude that the Bethe Ansatz is consistent only if the particles are indistinguishable, i.e., 
Aj is chosen such that either Ti = I or Ti = —I, and in these two cases we can compute all 
coefficients Ap from Aj as 

Ap = Wp{k)Ai (38) 

where yVp{k) is a product of operators Yi{kp^^-^ — kp.) and Z{2kp-^) in Eqs. (|36() and obtained 
by using repeatedly and . 

For Tj = +1 we get from Eq. (|36|) that Yi{u) = I independent of Ai, and we conclude that the 
momentum- dependent interaction is trivial for bosons. However, for Tj = — / (fermions) the Yi{u) 
are nontrivial phases. There are two different fermions cases, namely Ri = ±1. 

3.3 Duality 

It is interesting to note that there exists a simple duality relation between the fermionic pSp model 
considered here and the bosonic Cn delta- interaction model discussed in Section 2. Since the 
operators l^(u) and Z{u) for the latter model is identical with the ones of the fermions p6p model 
upon the substitution Ai — > 1/ci and A2 — > l/c2 (compare Eqs. (PT|) and for % = Ri = +1 
and Eqs. ^ and (EH) for % = Ri = -I), Eqs. ((23) and ^ imply that 

Ai.\ - ' — ap^p\ - - I'^q'i 
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where Ap are the coefficients of Section 2.2 and A^^ the ones in Section 3.2. This imphes that 
the bosonic wave functions of the delta model in Section 2.2 and the fermionic wave functions of 
the p6p-model in Section 3.2 are identical when restricted to the fundamental wedge 

A/: < xi < X2 < . . . < XAT, (40) 

provided that the coupling constants of these models are related as follows, 

1 1 
Ai = — and A2 = — . (41) 

Cl C2 

This can be seen also more directly: assuming that the eigenfunction ■0 of the Hamiltonian in Eq. 
is bosonic, Ti = R = I, it is enough to determine it in the fundamental wedge. Moreover, the 
continuity conditions in Eqs. (|llal b) are fulfilled automatically for boson wave functions, whereas 
the conditions on the derivatives simplify to 

{dxj - dn,.^^ - ci)V'U,=Xfc+o+ = 

(2a,^. - C2)V'U,=o+ = (42) 

for all X in the fundamental wedge. In a similar manner one finds that the fermionic eigenfunctions 
of the Hamiltonian in Eq. ©, Tj = ii = — /, are determined by the very same conditions in Eq. 
with ci^2 replaced by 1/Ai^2- 
This generalizes the duality previously observed in the A^-i case [4,5] to the Cn case. 



3.4 Physical interpretation 

As in the delta interaction case, one can restrict the eigenfunctions ip of the Hamiltonian in Eq. 
© by imposing the conditions in Eq. 1)28(1 . reducing the boundary conditions in Eqs. (|32al bl to 

{dx, -9^JV'U,=Xfc+0+ = {dx, - dxk)'4^\x,=Xk-0+ 

Mx,=Xk+o+ - V'U,=Xfc-o+ = 2Ai(aa;^. - jVL,=Xfe+o+ (43a) 

and 

9x, V'U,=o+ = for Rj = +1 ^^^^^ 

V'U,=o+ = 2A25a;^V'L,=o+ for Rj = -1 

where Xj > 0. This shows that the eigenfunctions of the Cn Hamiltonian in Eq. Q with the 
restriction in Eq. (|28() are identical to the ones of the ^at-i Hamiltonian 

N 

Ho = -Y. < + 2Ai Y.(^x, - d,,)6{xj - xu){d^^ - 9, J (44) 

j=l j<k 



restricted to the half-line, Xj > 0, and the boundary conditions at the origin given in Eq. (j4| 

Moreover, as shown in Appendix C.2, the eigenfunctions ifj above restricted to Xj > become 
identical to the ones of the Hamiltonian Hq + V{xj) on the full real line, Xj S M, but with an 
external potential 

V(r^ = / ^oe(-x) + ^odx5{x)d^ if R, = +1 . . 

I VoQ{-x) + 2\2dJ[x)d:, \iRj = -l ^ ' 

in the limit Vq — > 00. 
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4 Concluding remark 



As discussed in the Introduction, there exists a 4-parameter family of local interactions [3], and 
the delta- and p^p-interactions only correspond to one-parameter subfamilies each. It is therefore 
natural to ask: What about the other local interactions? Are there other cases leading to exactly 
solvable models? It is thus interesting to note that there is a simple physical interpretation of the 
four parameter family of local interactions which seems more natural than the ones given before [3]: 
in the simplest case they correspond to the following generalization of the Hamiltonians in Eqs. 
(HJ) and 0, 

H = -dl + c5{x) + Aa,5(x)5, + -id^5{x) - 75(x)5,, (46) 

which obviously is the most general hermitian Hamiltonian with interactions localized in a; = 
and containing only derivatives up to second order (higher derivatives than that do not lead to 
physically acceptable boundary conditions). This Hamiltonian is formally self-adjoint for arbitrary 
parameters c, A G M and 7 G C, and it indeed corresponds to the 4-parameter family of local 
interactions mentioned above [14]. All these models have natural generalizations to the many- 
body case, but there is only one case where these latter models are known to be exactly solvable 
even for distinguishable particles by the coordinate Bethe Ansatz: (c, A, 7) = (c, 0, 0). It would be 
interesting to know if there are other exactly solvable cases. We plan to come back to this question 
elsewhere [14]. We only mention here that the many-body generalization of the Hamiltonian in 
Eq. ()46() describes identical particles only if 7 = 0, and to find exactly solvable case for non-zero 
7 therefore requires an extension of Yang's method of solution [2] (which only works for identical 
particle models). 
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Appendix A. Verification of consistency relations 

In this appendix we sketch the verification of the consistency relations in Eqs. H26al b) (Yang- 
Baxter relations and the Reflection equation). 



A.l Delta-interaction 

We start by writing the operators Yi in the following way: 

Yi{u) = a{u) + h{u)fi (Al) 

where 

a{u) = b{u) = (A2) 

m — Cl m — Cl 



Inserting this expression into the equations in H26a|) results in a number of relations between the 
coefficients a{u) and 6(u), one for each equation and different permutation operator. Most of them 
are trivially fulfilled, but the following ones are non-trivial: 

a[—u)a{u) + b{—u)b{u) = 1 

a{-u)b{u) + b{-u)a{u) = (A3) 
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and 



h{v)a{u + v)a{u) + a{v)a{u + v)h{u) = a{u)h{u + v)a{v). (A4) 

Inserting a{u) and b{u) from Eq. HA2|) they can be verified by straightforward calculations. To 
verify Eq. H26b|) we write the operator Z as 



Z{u) = a{u) + h{u)Ri (A5) 

where 

a{u) = h{u) = (A6) 

m — C2 m — C2 

Substituting this and Eq. HA1|) leads to the following non-trivial relation, 

b{2v)b{u + v)d{2u)a{u - v) + b(2v)a{u + v)a{2u)b{u - v) + 

+ a{2v)a{u + v)b{2u)b{u - v) = a{u - v)b{2u)b{u + v)a{2v) (A7) 

in addition to 

a{—u)a{u) + b{—u)b{u) = 1 

a{-u)b{u) + b{-u)a{u) = 0, (A8) 

the validity of which follow from straightforward calculations. 

We conclude that the Bethe Ansatz is consistent even for distinguishable particles. 

A. 2 Local momentum-dependent interaction 

In this case we get Yiiu) as in Eq. HA1|) but with 

lu — 1/Ai m — Ai 

With that the two equations in (|A3() hold true but the equation in HA4() does not. We therefore 
conclude that the Bethe Ansatz is not consistent for distinguishable particles. 

For indistinguishable particles we have Tj = ±7 and the Yang-Baxter relations in Eq. (|26a() 
are trivially fulfilled. Moreover, in this case it is also easy to check that the relations H26b|) hold 
true for Ri = ±1. 

We conclude that the Bethe Ansatz is consistent in the indistinguishable particle case but not 
in general. 



Appendix B. Representations of the group Wn 

In this appendix we discuss the irreducible representations of the group Wn = (Z/2Z)^ ys Sn- In 
particular we will show the following. 

Fact: There exists a set of irreducible representations of Wn isomorphic to the irreducible rep- 
resentations x± Pj where x± is a character (irreducible representation) of the (normal) abelian 
subgroup (Z/2Z)^ such that x±{Rj) = =tl /o?^ a/Z j = 1, 2, . . . , (same sign for all j) and p is 
an arbitrary irreducible representation of the permutation group S]\f. 
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To show this we wih use the notion of induced representations, following Section 8.2 of [15]. 
We start by determining the group of characters X = Hom((Z/2Z)^, C) of the sub group (Z/2Z) . 
The fact that it is generated by the reflections Rj obeying the relations (see e.g. page 21 in [13]) 

R] = I, j = l,2,...,iV (Bl) 

implies that the characters x 6 ^ are functions such that 

x{Rj) = e'"^", G Z (B2) 

for all j = 1, 2, . . . , A^. The group Wn acts on these characters by 

{wx){R) =x(.w-^Rw), yw eWN,X^ X,Re {Z/2Z)^. (B3) 

We now determine the orbits of the action of Sn in X, represented by a set Xi where i S X/Sjy- 
Using the fact that the adjoint action of Sn permutes the reflections Rj, Tj^RjTjk = Rk with Tj^ 
the transposition interchanging j and k, we conclude that the orbits of Sn in X can be represented 
by the characters 

where j,k = 1,2, . . . , N. For each i let {S]\f)i be that subgroup of Sn consisting of all P G Sn 
such that Pxi = Xi-, and let further Wi = (Z/2Z)^ • {SN)i- The structure of Xi implies that 
{SN)i = Si X Sn-i- The character Xi can be extended to all of Wi by setting 

Xi{RP) = xiR), R e (Z/2Z)^, P E {SNh (B5) 

Now let Pi be an irreducible representation of {SN)i and combine it with the canonical projection 
Wi {SN)i to yield an irreducible representation pi of Wi. By taking the tensor product of Xi 
and Pi we can now construct a set of irreducible representations Xi ^X" Pi of Wi. We denote the 
corresponding induced representation of the whole of Wn by Oi^p^. It follows from Proposition 25 
in [15] that all irreducible representations of Wn are isomorphic to such a representation Oi^p^. In 
particular setting i = and i = N we arrive at the claim stated in the Fact at the beginning of 
the section. 



Appendix C. Physical interpretation of boundary conditions 

In this appendix we substantiate the physical interpretation of the boundary conditions of the C n 
models given in Sections 2 and 3 in the main text. 

C.l Delta-interaction 

We first recall the eigenfunctions ip of the one particle Hamiltonian in Eq. (^. Since this Hamil- 
tonian is invariant under the reflection x —x these eigenfunctions can be chosen such that 
^l^{x) = ziz'tp{—x) = i^±{x), and they can be computed using the Ansatz 

^ / e-^'^ + ^±ei'^ forx>0 

'^^^^^ " 1 ± (e^^- + A±e-i^-) for X < ' 

and the boundary conditions in Eq. Q determine the constants A± as follows, 
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with A_ being independent of c corresponding to the fact that the delta interaction is trivial (i.e. 
invisible) for fermions. Obviously, these eigenfunctions obey 

-i>'l{x) = k'^ip+{x) for X > and ^^'(0+) = (c/2)V'(0+) (C3) 

and 

-tl;'^{x) = k'^'il^_{x) forx>Oand ^-(0+) = 0, (C4) 

which is the simplest non-trivial case = 1 of the general relation between the Cn model and 
the Ajy-i model discussed in Section 4.1. 

We now show that these eigenfunctions ip±{x) for x > are identical to the ones of the 
Hamiltonians 

H± = -dl + V^Q{-x) + g^b(x) (C5) 

with 

3+ = c/2 - and 5- = (C6) 

in the limit Vq ^ cxd. To show this we determine the eigenfunctions <^± of }i± with the Ansatz 

f e"'*^^ + S±e^'^^', for X > , , 

= \ C7±e-^ for x < ' ^^^^ 

and by straightforward computations we find 

B^ = p:S^^±^ and a; = ^/v^ (C8) 

for Vq > k"^. We thus see that 

A± = lim B± (C9) 

provided that g± are chosen as in Eq. ()('6|) . This shows that the eigenfunctions 0+ of the Hamil- 
tonian Hj^ on the full line in the limit Vq ^ oo become equal to il^+{x) for x > (and zero 
otherwise), and similarly for (/>_, and 

This computation substantiates the physical interpretation of the Cn model in case = 1. 
However, since this interpretation only involves the boundary conditions at Xj = which are not 
affected by the inter-particle interactions, this argument immediately generalizes to the > 1 
particle case. 



C.2 Local momentum dependent interaction 

The discussion for the Hamiltonian in Eq. ((71) is completely analogous to the one for the Hamil- 
tonian in Eq. given above, and we therefore only write down the formulas which change. 

Eq. ()('H) determining the even and odd eigenfunctions il)± remains the same but Aj^ and A- 
are (essentially) interchanged, 

where now the boson eigenfunction is unaffected by the interaction. Moreover, these eigenfunctions 
solve the following problems on the half axis, 

-i^l{x) = k^ij+{x) forx>Oand ^^(0+) = (Cll) 

and 

-i;'L{x) =k^ij-{x) forx>Oand ^-(0+) = 2AV'-(0+). (C12) 
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The physical interpretation of these boundary conditions is provided by the following Hamiltonians 
with external fields, 

H± = -dl + Foe(-x) + ~9±d^b{x)d^ (C13) 

which has eigenfunctions as in Eq. (jnH) but with 

B± = TT ^ -. and uj = VVo- k^, C14 

IK — LU/{1 + ujg±) 

which converge to A± for Vq ^ oo provided that, for example, 

= y% and g- = 2A. (C15) 
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